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Abstract 

Well-known causality arguments show that events occurring during or at the end of 
inflation, associated with reheating or preheating, could contribute a blue component 
to the spectrum of primordial curvature perturbations, with the dependence k 3 . We 
explore the possibility that they could be observably large in CMB, LSS, and Lyman-a 
data. We find that a k 3 component with a cutoff at some maximum k can modestly 
improve the fits (A% 2 = 2.0, 5.4) of the low multipoles (£ ~ 10 — 50) or the second peak 
{I ~ 540) of the CMB angular spectrum when the three-year WMAP data are used. 
Moreover, the results from WMAP are consistent with the CBI, ACBAR, 2dFGRS, 
and SDSS data when they are included in the analysis. Including the SDSS galaxy 
clustering power spectrum, we find weak positive evidence for the k 3 component at the 
level of A% 2 = 2.4, with the caveat that the nonlinear evolution of the power spectrum 
may not be properly treated in the presence of the k 3 distortion. To investigate the 
high-/c regime, we use the Lyman-a forest data (LUQAS, Croft et at, and SDSS Lyman- 
a); here we find evidence at the level A% 2 = 3.8. Considering that there are two 
additional free parameters in the model, the above results do not give a strong evidence 
for features; however, they show that surprisingly large bumps are not ruled out. We 
give constraints on the ratio between the k 3 component and the nearly scale-invariant 
component, r% < 1.5, over the range of wave numbers 2.3 x 10 -3 Mpc -1 < k < 
8.2 Mpc -1 . We also discuss theoretical models which could lead to the k 3 effect, 
including ordinary hybrid inflation and double D-term inflation models. We show that 
the well-motivated k 3 component is also a good representative of the generic spikelike 
feature in the primordial perturbation power spectrum. 



1 Introduction 

Inflation has become a cornerstone of modern cosmology; it not only solves critical cosmo- 
logical problems, but also provides possibilities of exploring the infant universe [H El E]. A 
crucial aspect of inflationary theory is the primordial perturbation power spectrum, which 



1 



connects the quantum fluctuations in the early universe to the formation of structure at 
later times. Inflation generically predicts a nearly scale-invariant primordial power spec- 
trum, which agrees well with cosmic microwave background (CMB), large-scale structure 
(LSS), and Lyman-a forest observations, such as the Wilkinson Microwave Anisotropy Probe 
(WMAP) [HE], Cosmic Background Imager (CBI) [6], Arcminute Cosmology Bolometer Ar- 
ray Receiver (ACBAR) [7], Two-degree-Field Galaxy Redshift Survey (2dFGRS) 0[9], Sloan 
Digital Sky Survey (SDSS) USUI], and Lyman-a forest (Viel et al. [El EH [T4] and SDSS 
Lyman-a [T5l [16] )- 

While producing a nearly scale-invariant spectral index has become a major criterion of 
model selection, a deviation from pure scale invariance is consistent with the current data, 
and in fact fits better than n s — 1. Further elaborations have been investigated; a well- 
known example, suggested by the first-year WMAP data (WMAP1) [17J, is the running 
spectral index model, which has a large and negative running of the spectral index and a 
large tensor-to-scalar ratio (although the evidence for this is weakened when the three-year 
WMAP data [WMAP3] is used [5]). Detailed investigation revealed that a partially running 
spectral index model, which has a constant spectral index on large and small scales but a 
running spectral index on the relevant scales (about 3.4 e-foldings), provides as good a fit 
to WMAP1 as the full running spectral index model [18]. Because of the limited range of 
scales, this could be described as adding a localized feature to the primordial spectrum, a 
topic which has also received considerable attention [THl EDI Ell E2] , and which is the theme 
of the present work. 

Other hints of peculiarities in the CMB data have prompted the investigation of spectral 
features going beyond a simple power law. For example, the reconstruction of the primordial 
power spectrum directly from the CMB angular spectrum data indicates the prominent 
feature of an infrared cutoff on the horizon scale [231 121] ■ There are also a few outlying 
multipoles at smaller scales, I ~ 10 — 50 and I ~ 540, which have inspired modifications 
to the spectrum, such as would be provided by introducing a sharp step in the inflaton 
potential [T7J [191 EB (2TJ. This introduces oscillations into the spectrum which allow for 
better fits to the data [251 121]. More recent analyses of this subject can be found in Refs. 
[23I2HII2SIEDIEIII321E31E11S51ESIEZIES1. Such improvements of course come at the price 
of adding more parameters to the model; questions of significance can be handled through 
quantitative model selection criteria [3S1 SHI El 142] . 

Because the current data do not justify adding a large number of parameters to the 
description of the power spectrum, it is important to consider models with a small number 
of extra parameters, and preferably with a strong theoretical motivation. In the present 
work, we will investigate a new kind of spectral feature in this category, namely the addition 
of a component which scales like fc 3 , in contrast to the nearly scale-invariant spectrum ~ k°. 
As we will describe in Section [21 this behavior was a generic prediction [43] based on the 
requirement of causality, prior to the invention of inflation. Our observation is that extra 
contributions to the spectrum arising after inflation, for example during reheating, could be 
expected to scale like k 3 at low k. In fact, it was recently shown that this effect can arise in 
hybrid inflation models [HI [45] and double D-term inflation models [46J . 

Of course, other qualitatively similar distortions could arise in particular inflationary 
models, having k n spectra with n ^ 3. In this paper we focus on the n = 3 case because 
of the theoretical motivations mentioned above, and upon which we will elaborate further 
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below. However, for much of our analysis of the data, we would not expect greatly different 
results for nearby values like n = 2 or n = 4, for the reason that the large-A; behavior 
must be cutoff at some scale k c ; then the feature of interest is a localized spike (since the 
data do not allow for larger features) which would exceed the underlying scale-invariant 
contribution to Vn{k) only over a rather limited range of k. In that sense, the k 3 model with 
a cutoff can be considered as a rather generic model for spikelike features in the primordial 
power spectrum. Nevertheless, we will compare the k 3 spike to k n spikes with n = 1,2,4,5, 
which will clarify this point. Spikelike features in the primordial power spectrum have been 
previously investigated in the literature. For example, Ref. [22J finds a strong spike in the 
smooth hybrid inflation model; its location, however, is beyond the current CMB-LSS scales 
and hence is not observable. In any case, we take the point of view that the k 3 model (with a 
cutoff) is sufficiently well-motivated to provide an interesting test for a feature of this kind, 
and we have done a state-of-the-art comparison of this model of spectral features with the 
CMB, LSS, and Lyman-a data. 

To describe the k 3 effect, we modify the primordial curvature power spectrum, Vn{k) = 
(L/2n) 3 ATrk 3 (\n ]l \ 2 ) (see Appendix EJ) , to the form 



V n {k) = V^k)+V 3 (k) 

ris — l 




+ V 3 (k 



n.-l 




(l) 



where k Q is a pivot point which we take to be ko = 0.002 Mpc -1 throughout (to be consistent 
with the WMAP collaboration). The amplitude ratio of the k 3 component to the nearly 
scale-invariant component is 

We refer to the model whose power spectrum is described by Eq. as the k 3 model. It has 
just one free parameter relative to the power-law ACDM model: the ratio r 3 . However, in 
any realistic model, the k 3 behavior must be cut off at some maximum scale k c ; for example, 
in tachyonic preheating, where the tachyon curvature rr? a is negative, only the perturbations 



satisfying k < a 



are amplified. For simplicity we introduce a sharp cutoff: if k > k c 



then the k 3 term in Eq. (jTj) is set to 0. We refer to this as the k 3 model. It has one more 
parameter, k c . 

In this paper, we examine the extent to which the k 3 model is consistent with the data. 
We begin in Section [2] by reviewing the theoretical motivations for this effect. Section [3] 
continues with a brief description of the CMB, LSS, and Lyman-a data used in the analysis. 
We investigate the evidence from WMAP in Section HI and explore the k 3 component on 
small scales by the high-£ CMB, LSS, and Lyman-a data in Section [5j In Section [6] we 
explore the parameter space of hybrid inflation and double D-term inflation models which 
could give parameters suggested by the previous phenomenological analysis. We compare 
the A; 3 spike to more general k n spikes in Section [3 We give conclusions in Section [HI In 
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Appendix |X] we clarify the relation between the k 3 spectrum (in modern notation) and the 
causality prediction made in the early literature. 



2 Motivation for the k 3 component 

Prior to the idea of inflation, it was argued on the basis of causality that any fluctuations 
which are created independently of each other must have a correlation function that vanishes 
at least as fast as k 3 as k — > 0. This was first proven in the context of Newtonian gravity 
and then extended to general relativity in Ref . [13] . In Appendix |A] we recapitulate (and 
translate into modern notation) this argument. The essential observation is that a density 
perturbation 5p(x) / p which consists of contributions from causal processes originating at 
positions x a has the form 

W =EFo(x _ Xo ), (3) 

P a 

where causality demands that the monopole and dipole moments of F a vanish in a volume 
extending over distances greater than the causal horizon [13 08]. In Fourier space, this 
implies that F a ^ must vanish at least as fast as k 2 for small k. Therefore the correlation 
(|5pk/p| 2 ) falls like /c 4 at small k. Appendix lAl shows how this corresponds to a k 3 spectrum 
for the curvature perturbation. 

Inflationary perturbations evade this constraint by having been in causal contact with 
each other before being driven out of causal contact by inflation. But perturbations produced 
at the end of inflation, or between two consecutive stages of inflation, for example during 
a phase of tachyonic preheating, should obey the constraint. A A; 3 component could lead 
to interesting effects at short scales, such as the production of primordial black holes. It 
could conceivably also distort the spectrum of CMB fluctuations at longer scales, if it is 
appropriately cut off at the short scales. Of course, there is no reason to believe that the A; 3 
growth should continue to arbitrarily large k, since it is derived in the region k — > 0, and 
unrestricted growth at large k would lead to an unphysical UV divergence in the power. The 
value of the cutoff k c beyond which the effect vanishes depends on the particular model giving 
rise to it. In the case of tachyonic preheating, the cutoff is determined by the mass parameter 
(m 2 ) of the tachyonic field, since only modes with physical wave number k 2 < a 2 \m 2 \ at the 
time of preheating undergo exponential growth. 

Explicit examples of k 3 spectral components have been discussed in Refs. [321 EDJ I5T] : 
they were generated by preheating at the end of chaotic inflation. It was shown that those 
models do not have enough freedom to give an observably large k 3 component once the 
COBE normalization is imposed. On the other hand, hybrid inflation models have more free 
parameters and thus have more likelihood to produce an observably large effect. In Ref. [IB] 
it was noted that double (D-term) hybrid inflation models E31 E3] can produce a spike 
with a k 3 spectrum and a cutoff, due to a tachyonic instability which triggers the transition 
between the two stages of inflation. More recently it has been shown [4"4"l |4"5] that tachyonic 
preheating even in the simplest model of hybrid inflation can generate an observably large k 3 
component or alternatively large non-Gaussianity (see also Refs. [55j [56j E3 EE] ) for certain 
ranges of parameters. This effect arises at second order in the cosmological perturbation, 
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and its importance is controlled by how large the fluctuations of the tachyon field can become 
before their back-reaction ends inflation. 

In Section [6] we will make a detailed study of the predictions of tachyonic preheating in 
the hybrid inflation for the k 3 perturbation and how it compares to the CMB and LSS data. 
In the next section we describe these data. 



3 Data 

To explore the experimental evidence for the k 3 component, we use CosmoMC, a publicly 
available Markov-chain Monte-Carlo (MCMC) engine for exploring cosmological parameter 
space [59] Q We use the first year and three-year WMAP data (WMAP1 and WMAP3) 
[H E] as primary data sets, analyzing them with the July 2005 and May 2006 versions of 
CosmoMC, respectively. For completeness and as a consistency check, we consider both 
WMAP1 and WMAP3 in this work. 

We include the high-^ CMB data, CBI [6j and ACBAR [7J, to explore the evidence for 
the k 3 component on small scales. The ranges of data we use are 300 < I < 3500 (14 data 
points) for CBI and 300 < £ < 3000 (13 data points) for ACBARl The first bands of both 
data are not included in the analysis since they are not well constrained. 

The LSS data, 2dFGRS [HI [9] and SDSS main galaxy sample [TO] JH are also used in our 
analysis. To use the linear theory, we set /c max = 0.15/i Mpc -1 in the galaxy-galaxy power 
spectrum as recommended by 2dFGRS and SDSS; their ranges are 0.022 Mpc" 1 < k/h < 
0.147 Mpc" 1 for 2dFGRS (32 data points) and 0.016 Mpc" 1 < k/h < 0.154 Mpc" 1 for SDSS 
(17 data points). SDSS also probes the nonlinear regime, which allows us to explore the 
possibility of fitting the nonlinear regime with the k 3 component. The ranges are extended 
to k/h = 0.205 Mpc" 1 (2 more data points). Further discussions can be found in Section 15721 
We also include the SDSS Luminous Red Galaxies (LRG) power spectrum [11]; the range is 
0.012 Mpc -1 < k/h < 0.087 Mpc -1 (14 data points) for linear theory, and it is extended to 
k/h = 0.203 Mpc^ 1 (6 more data points) for nonlinear theory. The results of 2003 [10J and 
2006 SDSS [11] data are referred to as SDSS1 and SDSS4, respectively, since they used the 
SDSS data release 1 and 4, respectively. 

Lyman-a forest data allow us to explore the high-/c regime. We use two different Lyman- 
a data sets to test the k 3 component. Viel et al. [T2] used the LUQAS sample [T3], which 
has z = 2.125, and its range is 0.0034 (km/s) -1 < k < 0.027 (km/s) -1 ; they also reanalyzed 
the result of Croft et al. [H] (z = 2.72) in the same range. The SDSS Lyman-a data [T5| [T6] 
have 2 < z < 4 and the range is 0.0013 (km/s)- 1 < k < 0.02 (km/s^EI The data sets we 



1 See http: / /www.cosmologist.info/cosmomc 



2 The offset lognormal matrix for 2000 + 2001 CBI was incorrectly incorporated into May 2006 
or earlier versions of CosmoMC [60] , so we use our offset lognormal matrix in this paper. See 
Ref. [60] for a discussion on the impact of using the wrong matrix; also see the CBI website, 
http://www.astro.caltech.edu/~tjp/CBI/data2004, for a link to the discussion on this issue. 



3 2002 and 2005 2dFGRS data are used along with WMAP1 and WMAP3, respectively; the differences of 
the results, however, are not significant. 

4 For the SDSS Lyman-a data, we do not use the default code in CosmoMC; instead, we use the patch 
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use, from CMB, LSS, to Lyman-a forest, cover k space from 10~ 4 Mpc -1 to 3 Mpc -1 and 
provide full sensitivity to a possible k 3 component on observable scales. 

In our comparison of the power-law ACDM model with the k 3 and k 3 models, we do not 
include the BB polarization power spectrum or the Sunyaev-Zel'dovich (SZ) effect. Most 
settings in CosmoMC are by default. Typical MCMC chains have on the order of 10 5 points 
and we search these chains to find the best fit points. 

4 Evidence from WMAP: Large Scales 

In this section we focus on the large scale (small k) region of the spectrum, using primarily 
the WMAP data to fit the k 3 model; for comparison we will also show the effect of combining 
with other CMB (CBI and ACBAR), LSS (2dFGRS or SDSS), or Lyman-a data. To get a 
significant effect at low k, the cutoff k c is necessary in order to avoid being dominated by 
the large amount of extra power at high k. Table [1] lists the best fit k 3 models for WMAP 
alone, and in combination with other CMB (CBI and ACBAR), LSS (2dFGRS or SDSS), or 
Lyman-a data. 

As can be seen in Table (TJ the k 3 model improves the fit to the data, relative to that of 
the power-law ACDM model, by a reduction in the x 2 of Ax 2 = 3.6 for the first year WMAP 
data, where 



Recall that Ax 2 = 4.7 for the running spectral index + tensor model [61 J, and so the k 3 
model is not as good as the running spectral index + tensor model when the first year 
WMAP data are used (both models have eight parameters). Nevertheless, the k 3 model 
gives a large amplitude ratio (r 3 (fc c ) = 0.69), and hence WMAP1 allows for a relatively large 
extra component in addition to the nearly scale-invariant spectrum, as shown in the top-left 
panel of Fig. [TJ More intriguingly, for the three-year WMAP data, we find an improvement 
of Ax 2 = 5.4; recall that A% 2 = 3.1 for the running spectral index + tensor model [61 J , 




and so the k 3 model gives a better fit than the running spectral index + tensor model using 
WMAP3. The top-right panel of Fig. [TJ shows the spectra; it is interesting that a small 
modification (rs(k c ) = 0.14) to the nearly scale- invariant spectrum can give such a large 



To better understand why the k 3 term improves the fits, the Cj T spectra are plotted. 
The lower part of Fig. [TJ shows the best fit Cj T spectra for the power-law ACDM model 
(solid, black) and the k 3 model (dashed, red) using the WMAP1 (left) and WMAP3 (right) 
data (dotted, green); the Cj E spectra are hardly changed, so they are not shown. It can be 
seen that both spectra are consistent over a large range, but the Cj T spectrum of the k 3 
model provides a better fit at low multipoles {I ~ 10 — 50) when WMAP1 is used. Therefore, 
the k 3 model can help explain the "glitch" in the low multipoles of the Cj T spectrum. It has 
also been argued that an exponentially increasing step [231 121] or oscillations [23 I2S] in the 
primordial power spectrum fit the low Cj T multipoles better than the nearly scale-invariant 
spectrum. Since the peak of the k 3 model is similar to an exponentially increasing step or 



provided by Anze Slosar. See |http: / / www.slosar.com/ aslosar/lya.h tml 



A X 2 = X 2 (ACDM)- X 2 (A; C 3 ). 



(4) 




A X 2 . 
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Table 1: Best fit models of fitting WMAP alone, and in combination with other CMB 
(CBI and ACBAR), LSS (2dFGRS or SDSS), or Lyman-a data. X 2 ' is the contribution to 
X 2 from all data excluding WMAP. 



Data 


x 2 /x 2 ' 

(ACDM) 


ArVAx 2 ' 


n s 
(ko) 


r 3 (k ) 




k c 

(Mpc- 1 ) 


^3 

(*c) 


WMAP1 


1428.8/- 


3.6/- 


1.005 


0.0911 




0.00394 


0.694 


+ CMB 


1463.1/30.9 


3.2/-2.5 


0.987 


0.0782 




0.00385 


0.561 


+ 2dFGRS 


1463.2/34.4 


1.7/0.0 


1.010 


0.0280 




0.00408 


0.236 


+ SDSS1 (linear) 


1445.9/16.9 


0.5/-0.2 


0.987 


0.0715 




0.00356 


0.407 


+ SDSS4 (linear) 


1443.3/12.5 


3.8/-0.4 


0.983 


0.0123 




0.00364 


0.747 


WMAP3 


11252.3/- 


5.4/- 


0.944 


1.62 x 10" 


-5 


0.0388 


0.139 


+ CMB 


11285.9/33.3 


4.1/0.4 


0.949 


2.38 x 10" 


-5 


0.0373 


0.178 


+ 2dFGRS 


11290.8/38.4 


3.7/-0.2 


0.947 


1.08 x 10- 


-5 


0.0383 


0.0890 


+ SDSS1 (linear) 


11274.0/18.9 


4.5/-0.1 


0.952 


1.26 x 10- 


-5 


0.0396 


0.113 


+ SDSS4 (linear) 


11264.9/12.2 


1.1/-0.8 


0.955 


1.12 x 10- 


-5 


0.0347 


0.0666 


WMAP1 

SDSS1 + 

WMAP3 


1452.5/23.6 
11278.1/24.3 


0.9/2.2 
0.9/2.2 


0.959 
0.943 


8.36 x 10- 
6.18 x 10- 


-7 
-7 


0.147 
0.131 


0.397 
0.222 


WMAP1 

SDSS4 + 

WMAP3 


1453.7/24.3 
11276.9/24.2 


1.6/4.1 
0.4/2.4 


0.950 
0.945 


24.0 x 10- 
9.76 x 10- 


-7 
-7 


0.0850 
0.0725 


0.223 
0.0567 


WMAP1 

Lj-a + 

WMAP3 


1454.1/24.9 
11279.9/26.9 


2.5/3.3 
2.5/3.8 


0.973 
0.959 


1.05 x 10- 
0.612 x 10 


-9 
-9 


1.49 
1.49 


0.521 
0.331 
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Figure 1: The best fit primordial power spectra (upper) and Cj T spectra (lower) for the 
power-law ACDM model (solid, black) and the k 3 model (dashed, red) using the WMAP1 
(left) and WMAP3 (right) data (dotted, green). The scales of the wave numbers k (Mpc -1 ) 
are shifted to match those of the multipole moment £. 



oscillations to some extent, it is not surprising that the k 3 model can offer an alternative 
physical explanation for the improved fits found in Refs. [23] EH ES, EH] . 

As for WMAP3, Tabled] and Fig. [I] show that the best fit k 3 peak appears around 0.04 
Mpc -1 , close to the second peak {£ ~ 540) of the Cj T spectrum (see the right panels of Fig. 
[TJ). This can be understood due to the failure of the best fit ACDM model to match the data 
within lcr in this region. Due to improved sensitivity of WMAP3 to these higher multipoles, 
the "glitch" at £ ~ 540 is statistically more significant to WMAP3 than to WMAP1, where 
the £ ~ 10 — 50 glitch took precedence. This explains why the best fit power spectra of 
WMAP1 and WMAP3 are very different (see Fig. [T]). We emphasize that these results are 
consistent with each other. One could have the k 3 peak appear at 0.004 Mpc -1 for WMAP3, 
but the resultant A% 2 is smaller (A% 2 = 2.0); similarly, one could adjust the k 3 component 
to fit the high multipoles of WMAP1, but the Ax 2 = 1-2 is smaller than that of fitting the 
low multipoles. (Of course, it is possible for a single data set to favor, to some extent, the 
simultaneous appearance of the k 3 feature at different scales; one of the objective of this 
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paper is to constrain such an effect, which will be given in Section EH) Obviously, a more 
complicated primordial power spectrum could have a better fit. In particular, introducing an 
additional k 3 peak at 0.004 Mpc -1 while keeping the k 3 peak at 0.04 Mpc" 1 fits better than 
either peak by itself. References [25l [26] make a similar observation; by changing the location 
and amplitude of the oscillations, their model is able to fit different ranges of multipoles. 
Unlike the oscillations which affect all multipoles beyond the second peak of the CMB power 
spectrum, the effect of the k 3 component for WMAP3 in the best fit model is localized at 
the second peak. 

Figure [2] shows the distributions of the parameters for the power-law ACDM model 
(black) and the k 3 model (red) using the WMAP3 data. Solid (ACDM) and dashed (k 3 ) 
lines are marginalized probabilities, and dotted (ACDM) and dot-dashed (k 3 ) lines are mean 
likelihoods of samples; for Gaussian distributions, they should be the same. The mean 
likelihoods indicate whether the parameters are really being constrained, or affected by the 
priors and the volume of samples; they also show how good a fit one can expect [59] . There 
are four Markov chains for the model, and they satisfy the Gelman and Rubin convergence 
test R — 1 < 0.03. The Gelman and Rubin "variance of chain means" / "mean of chain 
variances" convergence test generally demands that R — 1 < 0.1 for each parameter; for 
example, the WMAP collaboration demands that R — l < 0.1 [5]. However, smaller numbers 
usually indicate better convergence Jl 

One can see from Fig. [2] that introducing the k 3 component does not have significant 
statistical influence on the standard parameters. Nevertheless, there are inconsistencies 
between the marginalized probabilities and the mean likelihoods of samples at the higher 
ends of Q^h 2 , Q c h 2 , 6, n s , and lnP^. These inconsistencies, along with the small bumps of 
lnr3 and lnk c , indicate that there is another local best fit point in the parameter space, 
i.e. the k 3 model, which fits the low multipoles. This can be seen in the corresponding 
probabilities for WMAP1 (which we do not show): there the smaller bumps seen in the 
distributions for WMAP3 become the more significant features, and shift the most likely 
values. 

We also notice that the marginalized probabilities of In r 3 and In k c are far from Gaussian. 
This is related to the priors of lnr3 and lnk c : 

-25 < lnr 3 (fc ) < 0, (5) 
-8 < ln£; c [Mpc" 1 ] < -0.25. (6) 

To justify these priors, we note that the amplitude ratio lnrs(fco) = —25 gives a negligible k 3 
component even on small scales, i.e. r 3 (l Mpc -1 ) ~ 10~ 3 . On the other hand, if lnrs(fco) > 0, 
then the k 3 component will dominate over the primordial power spectrum and it should have 
been observed. The prior of lnfc c is set to cover observable k space: Ink [Mpc -1 ] = —8 
roughly corresponds to the largest observable scale and Ink [Mpc -1 ] = —0.25 is the highest 
value found to be called by CosmoMC when the CMB data are used. These are not true priors 
in any rigorous sense. From a phenomenological perspective these models are unbounded 
from below and any proper Bayesian estimate of parameters is likely to be swamped by the 
large part of the allowed parameter space that provides effectively no observed deviations 



See CosmoMC's website, http://www.cosmologist.info/cosmomc for further discussions 
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Figure 2: The distributions of the parameters for the power-law ACDM model (black) and 
the k 3 model (red) using the WMAP3 data. Solid (ACDM) and dashed (k 3 ) lines are 
marginalized probabilities; dotted (ACDM) and dot-dashed (k 3 ) lines are mean likelihoods 
of samples. (k = 0.002 Mpc^ 1 .) 



from scale invariance. This can be seen in the marginalized probabilities of Fig. [2j However, 
for understanding to what extent the data allow such deviations these priors permit practical 
exploration of the allowed parameter space. 

The choice of priors still allows for degeneracies. For any value of In r 3 in the prior, taking 
the minimum value of In k c makes the k 3 contributions negligible, and vice versa. This fact 
explains why the higher tails of the In r 3 and In k c distributions in Fig. [2] do not decay, since 
the k 3 model becomes the power-law ACDM model in that region and has a large volume 
of samples. Similarly, the lower ends of the In r 3 and In k c distributions have a larger range 
where the k 3 component is negligible and hence give rise to higher values of marginalized 
probabilities, as seen in Fig. [2j 

Plotting the two-dimensional marginalized surfaces makes the argument clearer. Figure [3] 
shows the two-dimensional marginalized probabilities (and their 68% and 95% CL contours) 
and mean likelihoods of samples of In r 3 and In k c for the k 3 model using the WMAP3 alone 
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(left), and in combination with SDSS4 (linear) and Viel et al. Lyman-a data (right). (See 
Section [531 for further discussions.) The waterfall shape surfaces demonstrate that the high 
In r% and In k c regions are ruled out by the data. The high plateaus are the regions having a 
negligible k 3 component and a large volume of samples. The peaks in the probabilities and 
mean likelihoods of samples represent the best fit models. We normalize the marginalized 
probabilities and mean likelihoods of samples at their maximum points and show them in 
logarithmic scales, so that one can see how far the plateaus are from the best fit points in 
terms of A% 2 ~ —2 lnL|§ As we discussed above, there are two peaks in the mean likelihood 
of samples, the smaller one corresponding to fitting the glitch at low multipolesJll Since the 
plateau is near the la confidence level, the projection of the two-dimensional marginalized 
probabilities shows many contours where the surface crosses la; see the bottom panels of 

Fig. m 

Table [2] lists the best fit models and marginalized values for the power-law ACDM model 
and the k 3 model using the WMAP3 data. Due to the points just discussed, the one- 
dimensional marginalized values do not give a good estimate or reliable upper limits on In r% 
and In k c . (However, we will give constraints on the maximum amplitude ratio in Section l5\4l ) 
Nevertheless, the improvements of fitting the irregularities at the low multipoles (£ ~ 10 — 50) 
or the second peak (£ ~ 540) could give an additional contribution as large as 14% to 69% 
to the nearly scale-invariant spectrum (see Table [I]). This observation not only provides 
evidence for the k 3 model but also demonstrates the extent to which the data still leave open 
the possibility of large deviations from scale invariance in the power spectrum. 

Furthermore, Fig. [2] shows that the best fit points can have significant variations from 
the ACDM values, as can be seen from the mean likelihoods of samples. With our choice of 
priors these points contribute little to the marginalized likelihoods, but it is surprising that 
adding localized features can move physical parameters such as the angular size of the sound 
horizon by order "lcr." With future polarization data this should be less of an issue, but the 
robustness of physical parameters to the physics of inflation could be important for future 
dark energy studies that rely on CMB-derived physical parameters. 

5 Fitting Small Scales 

To explore the evidence for the k 3 component on small scales, we now focus on other CMB, 
LSS, and Lyman-a data. The results are discussed in detail in the following subsections. 

5.1 Fitting High-f CMB or LSS Data 

Including the CBI and ACBAR data allows us to test for the possible presence of the k 3 
component in the high-£ regime. The best fit models are listed in Table [TJ As can be been, 
the best fit power spectra of joint WMAP, CBI, and ACBAR are close to those of fitting 



6 Some of the points at the bottom of the waterfall are zero, so their logarithm will give negative infinities. 
To solve this problem, we plot 2 ln(L + e) instead, where we take e — 4 x 1CP 6 so that the minimum values 
are —25, roughly 5a from the best fit points. 

7 The peaks are not obvious, since they are plotted on a logarithmic scale. 
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Figure 3: The two-dimensional marginalized probabilities (top; 21nL), mean likelihoods of 
samples (middle; 21nL*), and marginalized probability contours (bottom; 68% and 95% 
CL) of lnr3 and lnk c for the k% model using the WMAP3 alone (left), and in combination 
with SDSS4 (linear) and Viel et al. Lyman-a data (right). The dashed (red) lines in the 
marginalized probability contours give the best bounds on the amplitude ratio over the 
observable wave numbers (see Section 15.4)) . The surface plots have an artificial floor of 
21nL = —25 inserted for plotting purposes (see footnote [6]) . 
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Table 2: The best fit power-law ACDM model and the k 3 model, and their marginalized 
(mean, 68% CL) values for the WMAP3 data. 
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a See Section [5T4l for the constraints. 



WMAP alone. This implies that the improvements to the fits come primarily from the 
WMAP data (£ ~ 10 - 50 or £ ~ 540), and hence the CBI and ACBAR data do not seem to 
give any additional evidence for a k 3 component. Although there is some excess power in the 
CBI and ACBAR data at small scales, I > 1000, this cannot be attributed to the k 3 effect 
because Silk damping (also known as diffusion damping) kills any sensitivity of the CMB to 
the primordial power spectrum in this region. Any perturbation to the exponential damping 
tail seen in the CMB angular power spectrum can only arise from astrophysical effects long 
after inflation. Thus, no evidence for a k 3 component at large k can be obtained just using 
CMB data. 

We therefore consider the evidence for a k 3 component at large k using LSS data. Table 
[1] shows that the locations and amplitudes of the k 3 peaks are close to those inferred from 
the WMAP data alone. We see that adding a k 3 peak at large k (< 0.15/i Mpc -1 for 
2dFGRS and SDSS1 [linear]; < O.lh Mpc -1 for SDSS4 [linear]) fails to improve the fit to 
the galaxy-galaxy power spectrum. 

To further quantify the significance of the WMAP data, we compute a quantity x 2 ', 
defined to be the contribution to the total x 2 coming from all data except WMAP. These 
values are shown in Table [TJ Since the k 3 peaks at ~ 0.004 Mpc -1 (WMAP1) are outside 
the ranges of the CBI, ACBAR, 2dFGRS, and SDSS (linear) data which we use (see Section 
[3]), these data are not directly sensitive to such a k 3 component. Indirectly, however, they 
have an effect, since the best fit k 3 model of WMAP1 favors a bluer tilt than do the CBI 
and ACBAR data. This explains why the k 3 peak at low k is somewhat disfavored by these 
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data: A% 2 ' = -2.5. On the other hand, for the k 3 peak at ~ 0.04 Mpc" 1 (WMAP3), 
the high-£ CMB and LSS data exhibit no preference, giving A% 2 ' of order 0(0.5). Thus, 
although the CMB and LSS data do not give positive evidence, they are still consistent with 
a k 3 component which introduces a perturbation on the nearly scale-invariant primordial 
spectrum whose magnitude is 7% to 18%. 



5.2 Fitting the Nonlinear Regime 

The SDSS data also probe the nonlinear regime, and we include data up to k/h = 0.2 Mpc -1 , 
as recommended by the SDSS team. As reference models, we consider the power-law ACDM 
model taking into account nonlinear evolution effects; we also fit the k 3 model using the 
nonlinear theory. The total x 2 ' s ) the x 2 's of the galaxy-galaxy power spectra, and the best 
fit k 3 models are shown in Tabled! As can be seen from the table, the results of WMAP1 or 
WMAP3 plus SDSS1 are consistent, and both k 3 models fit slightly better than the power-law 
ACDM model. 

Figure H] shows the best fit galaxy-galaxy power spectra of WMAP + SDSS1 for the 
power-law ACDM and k 3 models using nonlinear theory (fitting to k/h = 0.2 Mpc -1 ). The 
vertical shifts between the results of WMAP1 and WMAP3 are probably due to the shift 
of n s from 0.959 to 0.943. The bumps at k/h ~ 0.18 Mpc -1 are believed to be the results 
of the k 3 peaks in the same regime. Therefore, our results show that, except for the low 
[k ~ 0.004 Mpc -1 ) and high (k ~ 0.04 Mpc -1 ) multipoles, the k 3 component could also 
appear in the nonlinear regime (k ~ 0.14 Mpc -1 ). However, since it is not theoretically well- 
understood how a power spectrum with a k 3 component goes nonlinear, further investigation 
is needed before drawing any firm conclusions. 
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Figure 4: The best fit galaxy-galaxy power spectra of WMAP + SDSS1 for the power-law 
ACDM model and the k 3 model (going up to k/h = 0.2 Mpc -1 ). The SDSS1 data are also 
shown. 

Table [U also shows the results of fitting WMAP and SDSS4. The k 3 component can 
fit the data well (A^ 2 ' = 2.4 for WMAP3 + SDSS4). However, the bumps caused by the 
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k 3 component are not as obvious as in the SDSS1 samples, so the power spectra of SDSS4 
are not shown. The total Ax 2 's are smaller in the case of SDSS4 because the LRG sample 
favors a lower cutoff, k ~ 0.08 Mpc -1 , where the WMAP data are better fit by a nearly 
scale-invariant spectrum. 

We emphasize that when the nonlinear SDSS data are included, one can still adjust the 
k 3 component to appear in the WMAP regime. For example, the global best fit model for 
WMAP3 + SDSS1 has Ax 2 = 2.1 (though Ax 2 ' = -0.4), which is slightly better than for 
the best fit model where the k 3 component appears in the nonlinear regime (Ax 2 = 0.4 and 
Ax 2 ' = 2.4). 



5.3 Fitting the Lyman-a Data 



The Lyman-a forest probes the highest-fc regime, and it usually has a large effective redshift 
(z > 2). Therefore, even though it is probing smaller scales, the scales are still in the linear 
regime because it is looking at early times (smaller scales go nonlinear first). 

As mentioned in Section [3j we use two different Lyman-a data sets to test the k 3 compo- 
nent: Viel et al. rT2J[T3J[TI] and SDSS Lyman-a [T^llIS] data. The lower part of Tabled] lists 
the best fit k 3 models of WMAP + Viel et al. Lyman-a. Figure shows the corresponding 
best fit linear power spectrum at z = 2.125 and 2.72, respectively. The results of WMAP1 
and WMAP3 are consistent except for the vertical shifts, which are probably due to the 
change of the spectral index from 0.973 to 0.959. In this case, the improved fit (Ax 2 ' = 3.8 
for WMAP3 + Lyman-a) is due to accounting for a bump at ~ 0.02 (km/s) -1 in the observed 
power, which corresponds to a k 3 peak in the primordial power spectrum near ~ 1.5 Mpc -1 . 
The best fit k 3 model for WMAP3 + Viel et al. has matter density fl m = flb + ^c = 0.27, so 
the correspondence between the two different ways of specifying the wave number is given 
by 1 (km/s)" 1 ~ 97h Mpc" 1 for z = 2.125 and 1 (km/s)- 1 ~ lOAh Mpc" 1 for z = 2.72 [62]. 
As in the case of WMAP3 + SDSS1 (nonlinear), the global best fit model of WMAP3 + Viel 
et al. has the k 3 peak appear in the WMAP regime, which gives Ax 2 = 5.7 (Ax 2 ' = 0.2). 
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Figure 5: The best fit linear power spectra of WMAP + Viel et al. Lyman-a for the power- 
law ACDM model and the k 3 model. The LUQAS (z = 2.125) and Croft et al. (z = 2.72) 
data are also shown. 
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We also test the k 3 . model using the SDSS Lyman-a data. It turns out that there 
is no similar evidence for the k 3 component as in the case of WMAP + Viel et al, i.e. 
the maximum ratio r 3 (£; c ) is of order o(0.01), so the k 3 model is not distinguishable from 
the power-law ACDM model and hence the A% 2 is of order O(O.l). At first glance it is 
a bit surprising that adding two new parameters does not improve the fit substantially. 
However, the SDSS Lyman-a data are not given in a band-power description of Ph{k)', 
instead, the information is extracted in terms of the amplitude, k 3 P^(k, z)/2n 2 , its tilt, 
n c s = d\iaP^(k)/dliak, and the running, a e s = dn e ff/d\iak, at a pivot redshift z p = 3.0 and 
a pivot wavenumber k p = 0.009 (km/s) -1 [TBI [IS] . It is not obvious that such a description 
of the data can accommodate the extra freedom available in the k 3 model. 

5.4 Constraint on the Amplitude Ratio 

In the absence of significant positive evidence for a k 3 distortion, we now consider the 
experimental limit on the magnitude of the k 3 component in the primordial power spectrum. 
The right panels of Fig. [3] show the two-dimensional marginalized probabilities (and their 
68% and 95% CL contours) and mean likelihood of samples of In r 3 and In k c for the k 3 model 
using the WMAP3, SDSS1 (linear), and Viel et al. Lyman-a data. There are three peaks 
in the mean likelihood of samples, and they correspond the local best fits of WMAP (low 
and high multipoles) and Lyman-a. But only the peak of WMAP (high multipoles) in the 
marginalized probability can be seen, since it has a large volume of samples. 

Again, the waterfall-like surfaces indicate that the high lnr3 and \nk c region is ruled 
out by the data, but the one-dimensional marginalized constraints on In r 3 and In k c are 
meaningless, because of the strong mutual dependence of the two parameters. However, the 
constraint on r 3 (k c ) can be inferred from the two-dimensional marginalized contours. From 
Eq. ([1]), the maximum amplitude ratio is 

so we have 

lnr 3 (fc c ) = lnr 3 (fc ) + (4 - n s ) \nk c - (4 - n s ) \nk 

~ lnr 3 (fc ) + 31nA; c — 31nfc , (8) 

where we have used the fact that the current best fit model has n s ~ 1. Therefore, plotting 
a series of lines with approximately the same r 3 (k c ), 

lnr 3 (A;o) cx — 3 In k c , (9) 

and finding their intersection with the two-dimensional marginalized contour gives the best 
bound on the maximum amplitude ratio. From the lower part of Fig. [3] (the dashed, red 
lines), we find that 

r 3 < 0.7 (95% CL) for 1.8 x 10~ 3 Mpc' 1 < k < 0.3 Mpc" 1 (10) 
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when using WMAP3 alone, and 

r 3 < 1.5 (95% CL) for 2.3 x 1(T 3 Mpc" 1 < k < 8.2 Mpc" 1 (11) 

when using WMAP3 + SDSS4 (linear) + Viel et al. Lyman-a. 

It is possible that for a heavy tachyon, the cutoff on its perturbations lies below the 
smallest currently observable scales. Hence, we also investigate the model without a cutoff, 
i.e. the k 3 model. From the point of view of the algorithm, this is equivalent to putting 
the cutoff below the smallest scales in the data or, equivalently, above the largest k value 
called by CosmoMC. We tested the k 3 model with all the joint data used in the previous 
sections. All of the results give a negligible amplitude ratio, i.e. lnr 3 — > —25 (the lower prior), 
making the model indistinguishable from the power-law ACDM model. This gives rise to 
the following observational constraint on the k 3 model. The upper limit of the magnitude 
of the k 3 term can be read directly from the bottom-right panel of Fig. [3l taking the largest 
cutoff value (the upper prior is k = 15 Mpc" 1 when Lyman-a is used) in the contour gives 

r 3 (15 Mpc' 1 ) < 11 (95% CL). (12) 

6 Theoretical Models 
6.1 Hybrid Inflation 

In the previous sections we have shown that bumplike features in the power spectra of the 
CMB and LSS can be explained by the addition of a k 3 contamination of the nearly scale- 
invariant contribution. We noted in the introduction that preheating at the end of chaotic 
inflation predicts k 3 components which are too small to be observed on the CMB-LSS scales 
[49l l50l 151] , but recent work has found that observable effects can be generated for some 
ranges of parameters in hybrid inflation theories [SJ 05] , which are defined by the potential 

m 

V{<p,v) = ±(o*-J)\^<f + £tf?o*, (13) 

where (p is the inflaton and a is the tachyonic field. 

Inflation starts with large values of <p and ends after the field- dependent tachyon mass 
parameter becomes negative: 

ml = -\v 2 + g 2 ip 2 < 0. (14) 

For certain ranges of the parameters v, A, and g, the fluctuations in a can become so 
large that their contribution to the curvature perturbation, at second order in cosmological 
perturbation theory, exceeds the usual first-order contribution due to the inflaton (p. The 
growth of fluctuations gets cut off at some number of e-foldings after the onset of the 
tachyonic instability, when their energy density starts to exceed the false vacuum energy 
that drives inflation. At this point, inflation ends and there is no further evolution in the 
curvature perturbation on superhorizon scales. Even if iV* represents a very short amount 
of time during which the fluctuations can grow in magnitude, they grow exponentially fast, 
and so it is possible to achieve a large effect. 
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Without going into all the details, one can still appreciate why it is possible to obtain 
a spectrum of the form k 3 with a cutoff. In fact, the k 3 spectrum is completely generic for 
massive fields, i.e. with m 3> H, where H is the Hubble scale during inflation. This can 
be found by solving the Klein-Gordon equation for linearized fluctuations of the fields in de 
Sitter space; it is straightforward to show that ((5a) 2 ) ~ mr 1 / d 3 k for heavy fields, whereas 
((Sip) 2 ) ~ H 2 J d 3 k/k 3 for light (m <C H) fields. However, because of the redshifting of 
the physical wave number, a k 3 component in the spectrum is usually negligible compared 
to the scale-invariant contribution. Only if the k 3 component is amplified can it become 
observable. Thus, the other relevant effect is the amplification of the tachyonic fluctuations 
by the instability, which causes modes with k 2 < a 2 |m^| to grow as ex$\b\J\m 2 \ — (k/a) 2 }. 
We see that there is naturally a cutoff k c = a\m a \ on the modes which get amplified. Of 
course, this scale gets stretched by the Hubble expansion which occurs after inflation. At 
the present time, the cutoff scale is therefore given by 

k c =\m a & N % (15) 

Hi 

where H and Hi are respectively the Hubble rate today and during inflation, and N e is 
the number of e-foldings of inflation since the horizon crossing, N e ~ 60. (In the present 
analysis, we do not assume that N e ~ 60; rather, we determine N e according to the actual 
reheat temperature, assumed to be of the order A 1//4 t> since Xv 4 is the false vacuum energy 
density which drives inflation.) Notice that if \m a \ = H iy k c /H is just the ratio of scales 
which crossed the horizon at the beginning and at the end of inflation, as expected. In order 
to match the bumps in the WMAP data, we need k c /H to be on the order of 16 or 162. 
Due to the large factor e N ° in Eq. (TT5I) . this is difficult to achieve. However, one should 
recall that m a is field-dependent (Eq. ( TH"j) ). and it passes through zero shortly before the 
end of inflation. For some values of parameters, m a can still be small at the end of inflation; 
moreover, N e can also be relatively small if the reheat temperature is low. It is therefore not 
obvious whether small enough values of k c can be achieved to match the phenomenological 
requirements that were suggested by the previous sections. 

To investigate whether the parameters suggested by our fits can correspond to hybrid 
inflation, we have adapted the code used in Refs. [HJ US], which numerically computed the 
magnitude r 3 of the k 3 contribution, scanning over a large range of the model's parameters, 
g, A, and v. In the left panel of Fig. [6] we plot the values of r 3 (at k = 1 Mpc -1 ) versus k c 
obtained by scanning over the parameter values log 10 v = —1, —3, —5, —7, —9, log 10 g 2 > —30, 
and log 10 A < 0. We find no cases where the peak parameters k c and r 3 are in the range 
needed for the WMAP glitches at i ~ 10 — 50 and I ~ 540 (the diamonds in the figure). The 
left panel of the figure shows that points having small enough values of k c all have overly large 
values of r 3 . However, there are realizations in which k c is close to the 1 Mpc -1 suggested 
by the Lyman-a data, as well as points where k c is larger and is effectively removed from 
the analysis (the k 3 model). Such points are contained in the narrow box shown in the lower 
right-hand corner of the left panel of Fig. [6J 

The promising points in the left panel of Fig. [6] (those falling within the box) correspond 
to values of the hybrid inflation parameters shown in the right panel of the figure. These are 
shown as boxes in the g-X plane, for the values of log 10 i; = —1, —3, —5, —7, —9, along with 
excluded regions based on the previous analysis of Refs. [HJ ES] , which demanded that the 
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Figure 6: Left: values of r 3 (at k = 1 Mpc -1 ) and k c which are achieved by varying parameters 
of the hybrid inflation model. The diamonds are the best fit models of WMAP at low and 
high multipoles; the boxed area is the Lyman-a preferred region — 2 < log 10 r 3 < 1 and 
6 < log 10 (fc c /ifo) < 10. Right: preferred points (square) next to excluded regions (circles) 
for fitting Lyman-a data in the hybrid inflation model. The numbers (—1, —3, —5, —7, —9) 
indicate the value of log 10 v for the respective regions. 

k 3 effect not be too large. As expected, the points suggested by our current fits to the data 
are close to the boundaries of the excluded regions. 

6.2 Double D-Term Inflation 

We have seen that pure hybrid inflation is able to give a sizeable k 3 component, but with a 
cutoff k c larger than the ranges observable in LSS and Lyman-a data. A more complicated 
but still well-motivated class of models is able to bring k c into the observable part of the 
spectrum, as was shown by Ref. jl6]l§ These are double inflation models of the D-term 
hybrid type [521 [53], [M] , based on the superpotential 

W = a AA + A_ + f3BB + B^, (16) 

with six superfields A, A±, B, B±, with charges (cja-i 9b) under two U(l) gauge groups U(1)a 
and U(1)b, given by (0, 0) for A, B, (±1, 0) for A±, and (±1, ±1) for B±. It leads to a scalar 
potential of the form 

V = f (U ~ l\C\ 2 ) 2 + \pB*\Cf + 9 f (& - i|C| 2 ) 2 , (17) 

plus one-loop corrections (depending on A, B, C), where £j are Fayet-Iliopoulos terms, A and 
B are proportional to the modulus of the scalar component of the corresponding superfields, 
and C is the complex scalar component of B_. 

From the form of (IT7|) it can be seen that the field C becomes tachyonic when B rolls 
down to a critical value. However, in this model, the instability of C need not trigger the end 



'We thank Julien Lesgourgues for calling our attention to this work. 
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Figure 7: A sample power spectrum from double D-term inflation, adapted from Ref. |46j. 
The bold dashed line is the k 3 contribution from tachyonic fluctuations. The model param- 
eters were y/U = 3 x 10~ 3 Mpi, = 4.2 x 10" 3 M P i, g A = g B = lCT 2 , and (3 = 0.5 x 10~ 3 . 
For comparison, a k 3 model is also shown, with parameters r 3 (k ) = 1.3 x 10~ 3 and 
k c = 0.3h Mpc" 1 . 

of inflation; instead, there can be a second stage via rolling of the field A. The fluctuations 
due to the tachyonic preheating are found to have a k 3 spectrum, and a sharp cutoff at 
the scale which crossed the horizon during the phase transition. A k 3 spike is generated, 
as shown in Fig. UJ In this example, the spike is located at rather high k values, but there 
is plenty of freedom in the model to move this to arbitrarily lower scales. One merely 
needs to adjust the duration of the second stage of inflation so that the phase transition 
occurs sufficiently soon after horizon crossing of modes which are at the present horizon. 
This model is therefore a good candidate for generating the bumps in the low multipoles 
which we have investigated in Section |H We also show a k 3 power spectrum in Fig. [7J By 
appropriate choices of the parameters of the double D-term inflation, the wavenumber and 
the amplitude of the feature can be adjusted. However, the most closely corresponding spike 
in the k 3 model has a narrower width, so one would have to add another parameter to the 
k 3 ansatz so as to accurately mimic the prediction of the double D-term inflation scenario. 

7 k n Features in the Primordial Power Spectrum 

As mentioned in the previous sections, the k 3 model is theoretically well motivated, and 
we take the point of view that it is a good representative of the generic spikelike feature in 
the primordial power spectrum. To clarify this point, we compare the k 3 model to the k™ 
models with nearby values (n = 1,2,4,5). Since the evidence from LSS and Lyman-a data 
is moderate, we focus on the WMAP3 data. Table [3] shows the best fit models for different 
k n spike power spectra. 

It can be seen from the table that all fc™ models are consistent: a k n spike at k ~ 
0.004 Mpc -1 can fit the low-multipole glitches (£ ~ 10-50) and a k n spike at k ~ 0.04 Mpc" 1 
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Table 3: The best fit fc" models of WMAP3, fitting to the low multipoles {£ ~ 10 - 50, 
upper) or the second peak {£ ~ 540, lower) of the CMB angular spectrum. 



Model 


A X 2 


n s (k ) 


k c (10~ 2 Mpc" 1 ) 


r n (k c ) 


k 1 


0.2 


0.970 


0.434 


0.159 


k 2 


1.2 


0.985 


0.384 


0.446 


k 3 


2.0 


0.969 


0.376 


0.397 


k 4 


1.5 


0.974 


0.382 


0.271 


k 5 


1.9 


0.981 


0.369 


0.797 


k l 


5.8 


0.944 


3.94 


0.123 


k 2 


4.3 


0.942 


3.82 


0.153 


k 3 


5.4 


0.944 


3.88 


0.139 


k* 


4.8 


0.946 


3.80 


0.200 


k 5 


5.4 


0.959 


3.79 


0.191 



can improve the fits of the second peak (£ ~ 540) of the CMB angular spectrum. While all 
k n spikes have roughly the same order of improvement (Ax 2 ), it is interesting that the k 3 
model is slightly favored by the data. 

We emphasize that the power n is fixed at n = 3 in most of the investigation; only in 
this section so we allow it to change to some nearby values, to show that n = 3 is a good 
representative of the generic spikelike feature. 

8 Conclusions 

We have examined a class of perturbations to the nearly scale- invariant spectrum, coming 
from preheating, which is motivated by very general principles of causality. It gives a k 3 
component which adds to the usual k ns ~ x spectrum, up to some cutoff k c which is determined 
by the microphysics of preheating. We did a phenomenological analysis of this kind of 
distortion to the primordial power spectrum, using the Monte-Carlo Markov-chain algorithm 
provided by CosmoMC, and found that such a component could improve either the fits to 
the irregularities of the low multipoles (£ ~ 10 — 50) or the second peak {£ ~ 540) of the 
CMB spectrum, giving A% 2 = 3.6 or 1.2 at the low multipoles and A% 2 = 2.0 or 5.4 at the 
second peak, using the WMAP1 or WMAP3 data, respectively. Moreover, the amplitude of 
the k 3 contribution was found to be as large as 0.69 or 0.14 of the k ns ~ x part when fitting 
the low multipoles or the second peak, respectively. These results provide an intriguing 
suggestion for such a k 3 component, but more generally they indicate that the WMAP data 
are consistent with sizable deviations from a nearly scale-invariant spectrum. 
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We also studied the CBI and ACBAR data to investigate whether the k 3 component could 
explain excess power in the high multipoles seen by those experiments. Our results show 
that the CBI and ACBAR data are consistent with a k 3 component at the low multipoles or 
the second peak, but they do not give evidence for a k 3 component at the high multipoles; 
in retrospect this had to be the case, due to Silk damping. However, this limitation does not 
apply to LSS data, i.e. 2dFGRS and SDSS, which explore a similar range of k space as CBI 
and ACBAR. Again, the LSS data are consistent with the results from WMAP, but there is 
no evidence for the k 3 c model at k/h < 0.15 Mpc" 1 (0.1 Mpc" 1 for SDSS4). 

We also tested the k 3 model in the nonlinear regime of the galaxy-galaxy spectrum, k/h > 
0.15 Mpc" 1 (0.1 Mpc" 1 for SDSS4). We found that the k 3 component can tune the shape of 
the SDSS galaxy-galaxy power spectrum, giving an improvement of A% 2 = 0.9 (Ax 2 = 2.2) 
for SDSS1; SDSS4 favors a lower cutoff, k ~ 0.8 Mpc" 1 , which gives A X 2 = 0.4 (Ax 2 ' = 2.4). 
However, we feel that the implementation of nonlinear evolution in the likelihood code for 
SDSS has not yet been tested thoroughly in conjunction with a nonstandard spectrum such 
as the one we are using, so we reserve judgment as to these particular results. They should 
be taken as motivation for a more detailed study of the nonlinear regime. 

A further handle on the power spectrum at high k is provided by the Lyman-a forest data. 
We found that the Viel et al. data allow a large amplitude ratio (0.33), with an improvement 
of x 2 — 2.5 (Ax 2 ' = 3.8). The SDSS Lyman-a data, however, do not give similar evidence 
in the same regime. 

In brief, the CMB, LSS, and Lyman-a data are consistent with a nearly scale-invariant 
spectral index plus a k 3 component. We further showed that the k 3 spike is a good rep- 
resentative of the generic spikelike feature in the primordial power spectrum by comparing 
different k™ models. By adjusting both the location and the amplitude of the k 3 component, 
one can of course always find better fits. We emphasize that due to the addition of two 
free parameters, the evidence for the k 3 component is not compelling. However, our results 
show that the k 3 contamination is not ruled out by the data, even when its amplitude is sur- 
prisingly large. We have determined constraints on the magnitude of this extra component, 
finding an upper limit of r 3 < 1.5 (95% CL) on the amplitude ratio, over the range of wave 
numbers 2.3 x 10" 3 Mpc" 1 < k < 8.2 Mpc" 1 . 

We have also explored in some detail the parameter space of the hybrid inflation model, 
which was found to give an observably large k 3 component for some ranges of the model's 
parameters. The model is able to match the observations in the cases where the benefit of the 
k 3 component is less rigorously shown, namely in the highest-A; regions of the spectrum. For 
the features at lower k, the double D-term inflation model discussed in Section 16.21 appears 
to be ideally suited for generating spikes at the large scales we have investigated here. 
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A Relation of the k 3 Spectrum to the Causality Con- 
straint 



Due to changes of notation over the years, readers of the original causality argument due 
to Abbott and Traschen [43] may not immediately perceive that their result implies the k 3 
spectrum for causal perturbations. In this appendix we clarify the relation. Reference [4"3] 
assumes that 5p(x)/p = J2a c aF a {^ — x a ), where for convenience c a is a random variable 
which obeys (c a Cb) = c 2 5 a b. The correlation function in position space is 



'^\ = c2Effl(x _, )fn(y _ Xfl) . 

, P P a 



In Fourier space, 



' £" * £" 
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k 2 , which is the essential 



where we assumed that each F a has Fourier transform F ak ~ |k| 
restriction due to causality. If there are enough random centers x a , then the sum will 
approximately give a delta function, 
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where L is a box size. Therefore, one has 
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The k A behavior is potentially confusing, but we will now show that this translates into k 3 
behavior for the power spectrum of the curvature perturbation 1Z. 
The relation between TZ^ and Sp k /p is [62] 
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where Di(Q m )/Q m is the growth factor for the matter perturbation and T(k, Q m o) is the 
transfer function. Since the growth factor does not have k dependence and the causality 
argument applies to the limit k — > where T(k) is roughly a constant, 
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On the other hand, the power spectrum is defined through 

V n (k)= ^\nk 3 (\n k \ 



(23) 



(24) 
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Equations (j23J) and fl2U) imply that 



8 h 



P 



2tt\ 3 fc ^ .,. 



Making a comparison with Eq. (j2ip . we see that 

Vn(k) ~ k 3 



(25) 



(26) 



if <5p has the assumed behavior. 

The typical onset value for T(k) to decrease is of order 0.01 Mpc -1 . Therefore, if the 
causality constraint is going to give a k 3 spectrum, then the k 3 component appears only at 
k < 0.01 Mpc -1 , otherwise the causality constraint will give a power greater than 3. Of 
course, the decrease of the k 3 component could be more complicated than a sharp cutoff. 
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